Matter— wave emission in optical lattices: Single particle and collective effects 
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We introduce a simple set-up corresponding to the matter-wave analogue of impurity atoms 
embedded in an infinite photonic crystal and interacting with the radiation field. Atoms in a given 
internal level are trapped in an optical lattice, and play the role of the impurities. Atoms in an 
untrapped level play the role of the radiation field. The interaction is mediated by means of lasers 
that couple those levels. By tuning the lasers parameters, it is possible to drive the system through 
different regimes, and observe phenomena like matter wave superradiance, non-Markovian atom 
emission, and the appearance of bound atomic states. 



Recent progress in atomic physics has allowed exper- 
imentalists to trap atoms in optical potentials at very 
low temperatures. This has led to the observation of 
several interesting phenomena in which atom-atom in- 
teractions play a predominant role. With atoms loaded 
in optical lattices it is nowadays possible, for example, 
to reach the strong correlation regime where quantum 
phase transitions between superfluid and insulator phases 
[l|, y , Tonks-Girardeau gases [3| , or even entanglement 
between neighboring atoms can be observed. Those ex- 
periments have triggered a large amount of theoretical 
work proposing and analyzing new experiments where 
intriguing condensed matter behavior could be observed. 
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In this work we show that with the same systems it 
is possible to observe a broad spectrum of different phe- 
nomena usually connected to light-matter interactions 
(see also [J, la, |6| for related setups). In our setup, the 
role of matter is played by the absence/presence of one 
atom in the ground state of an optical potential, whereas 
the role of light is played by weakly-interacting atoms in 
a different internal state which is not affected by the op- 
tical potential. The coupling between those two systems 
is induced by Raman lasers, which connect the two in- 
ternal states of each atom (see Fig. [T]). As we will show, 
the Hamiltonian that describes this situation is very sim- 
ilar to that describing the interaction between two-level 
atoms and the electromagnetic field within a photonic 
crystal (PC). By changing the laser and optical trapping 
parameters it is possible to drive the system to different 
regimes where a rich variety of phenomena can be ob- 
served. These include the spontaneous polarization of the 
system predicted by the mean field theory 7], collective 
effects in the emission of atoms from the lattice [8, 9? ], 
and the formation of a bound trapped-untrapped atom 
state, analogous to the atom-photon bound state that 
appear when atoms within a pho tonic crystal emit pho- 
tons within the gap region 10|) lUl lid, ll3l . Moreover, it is 
possible to reach a regime in which weakly confined atoms 
drive atom-atom interactions between strongly confined 
ones, giving rise to effective Coulomb-like interactions 
between them. 

We consider N cold atoms with a ground state hy- 



FIG. 1: (Color online) Schematic diagram: Lasers produce 
two-photon Raman transitions to an untrapped state. Fig 
A: For laser detunings A > 0, transitions are into a band of 
non-trapped states. Fig B: A < 0, transitions are into the gap 
region, and a trapped-untrapped atom bound state is formed. 



perfine level \a) and frequency w^ that is trapped by an 
optical lattice with M sites and lattice period do. The 
motion of the atoms is restricted to the lowest Bloch band 
in the coUisional blockade regime, where either one or no 
atom occupy each potential well of the lattice, which we 
will approximate by a harmonic oscillator of frequency 
ujo- Then, we can replace the creation operator at each 
site by a^ = |l)j(0|, which describes transitions from the 
Fock state |0)j with no atoms at site j, to a state |l)j 
with one atom f4|. The atoms have an additional inter- 
nal level, 1 6), that is not affected by the lattice potential 
(what can be achieved with state dependent potentials) 
and has a frequency ujh. We introduce the field operator 
'0b(r) = (l/VV) X^k^k^'*"^' where V is the quantization 
volume, and b^^ is the creation operator of an atom in \b) 
with momentum k. 



In a similar way as in an atom laser setup IJ], two 
lasers are then used to induce two-photon Raman tran- 
sitions between | a) and | b) . The lasers have a two-photon 
Rabi frequency fi, and their frequencies and momentum 
differences are ujl — uJi — loi and k/, = ki — k2 respec- 
tively. When tuning them close to a two photon reso- 
nance and far from single photon resonances, an effective 
Hamiltonian is obtained which in the interaction picture 



can be written as (S = 1) 

J k 

Here rj denotes the positions in the lattice, and Ak — 
k'^ /2m — A, with A = lol — {(^t — ^a) the laser detuning 
(and uja — uj^ + ojq/2). The coupling constants are gu = 
ne-^ofeV2(8^3/2;s^3/y)i/2^ where Xo = (l/2mu;o)i/" is 
the size of the wave function at each site. 

The similarity of Hamiltonian llj with that describing 
the interaction of atoms with the electromagnetic field 
is apparent. The dispersion relation of the atomic bath 
is contained in A^, and resembles that of the radiation 
field in a three dimensional and infinite PC near the 
band-edge [71, [ill, ll3[. Furthermore, in our set-up one 
can easily control several external parameters: fJ, which 
determines the coupling strength. A, which determines 
the resonance conditions, the number of atoms N and of 
sites M, the lasers wave vectors k^, and the dimension of 
the trap and the lattice. Thus, we expect to observe a 
rich variety of phenomena with our system, ranging from 
some well-known from the field of Quantum Optics, to 
other which are difficult to access in that field. We will 
start out with a single excitation [N = M = 1), and 
then consider collective effects (M > 1), for the different 
regimes dictated by the control parameters Q and A. 

An atom within a single trap constitutes the sim- 
plest setup but it still gives a very good insight into 
the problem. The wave function of the system has 
the form |*(i)) = A(t)|l, {0}> + Ek ^kW|0, Ik), where 
|1,{0}) describes the atom in the trapped state and 
no free atom present, and |0,lk) represents no atom 
in the trapped state and a single untrapped atom in 
the mode k. Using the Schrodinger equation we have 
A(t) = - /q dTG{t - t)A{t), where 



GW^Eff'^'^^"* 



„'i(At— arctan [a;oi]) 

^'- rT2 — ' (2) 



with Ut = 2^/1 -|- iuot, is the correlation function of the 
environment. An analytical solution can be obtained by 
assuming that trapped atoms are strongly confined so 
that LUo ^ ^, A. For a 3D bosonic field, this leads to 
Goo(i) = -ae'('^*+''/''V^^^^, which is singular at the ori- 
gin, but describes correctly times t 3> l/i^o- Except for 
the value of a = n?/ujQ , Gooit) is identical to the corre- 
lation function of the radiation field within an anisotropic 
PC, as described in [12| . In the same way, a ID environ- 
ment (produced by trapping the atoms in \b) in a ID 
harmonic trap), gives rise to a correlation function simi- 
lar to the one appearing for the radiation field in isotropic 
PCs. Using the Laplace transform method, we get 
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FIG. 2: Evolution of the atomic population |^(i)p in loga- 
rithmic scale for different detunings. Solid, dotted, dashed, 
dot-dashed and dot-dot-dashed lines correspond respectively 



to A/a-' = 



-1, —0.2, 0.2, 8 {ujQ = oo), where one can rec- 



ognize the regimes explained in the text. Inset [131: Steady 
state population |A^*p — |yl(co)|^. Solid, dashed and dotted 
lines correspond respectively to fl/ioo = 0.05, 0.025, 0.01. 



with I{a, A, t) = (ae-/V^) J^ dxj^^ 
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ing r± = -(a/2) ± ^(q/2)2 - A we have: (i) If a^ /2 > 
A > 0, then a = 0. (u) If A > a^/2, then n = r_ and 
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Depending on the parameters, we can have very differ- 
ent behaviors. For A > there is no trapped atom left 
in steady state, whereas for A < this is not the case 
(Fig. [2]). Thus, there is a quantum phase transition at 
A = analog to that found in the spin-boson model [l5| . 
Moreover, for A < 0, the emitted atomic field is in the 
form of evanescent modes exponentially localized around 
the trapped atom, what leads to the trapped atom- 
untrapped atom bound state. This follows from the prob- 
ability of finding a radiated particle at position r at a long 
time t, |*b(r,t)p_^_(cif^mX^/27rr)2e-i'"['=ol'-e-i'»['-?l*, 
where k^ = •\/2m(A — r^) is imaginary if A < 0. For 
A > we have two different regimes: for A ^ a^ we just 
have an exponential rate Fg = 2^ ^/t^KJus^^ whereas 
in the opposite limit the evolution does not follow such 
a law. The first corresponds to the Markovian regime, 
where the correlation time t^ — A"^ of the environment 
(here untrapped atoms) is shorter than the typical evo- 
lution time of the trapped atom (Fq ^) |l7|. 

We now study the dynamics of atoms in a lattice with 
M sites, choosing luq ^ ^, A. Guided by the previous 
analysis, we will consider the regimes A > (Markovian 
and non-Markovian) , as well as A < 0. 

In the limit where FcoUTc "C 1, where FcoU gives the 
typical evolution time of the trapped atoms, we can ana- 
lyze the problem under the Born-Markov approximation. 
The dynamics of the atoms in the lattice is dictated by 
the quantities 



A(i)=cie^('''+^)*-l-/(a,A,t), 



(3) 



|i-j| 



p-J'li-Jl/e-l-in-j-ki 

drGi_j(r)=z|Fo|e r^^. , (4) 



for i ^ j. Here, v = z, 1 for A > and A < 
0, respectively, and the correlatfon function G'i_j(t) = 



Gi_j(i) = G{t)e'^U/(^x>'t) 



(5) 



with G{t) given by ^. Similar to the radiative case, 
the coefficients r|i_j| describe the dipolar interactions be- 
tween the sites i and j. The quantity ^ = l/(|fcoMo)i with 
fco = v2mA J.16], quantifies the range of the interactions 
which, according to (j4|), has a Yukawa form. 

For A > 0, the situation under study resembles that of 
a set of M atoms in a lattice of constant do interacting 
with the electromagnetic field with resonant wave-vector 
ko, and where N corresponds to the number of excited 
atoms. Thus, phenomena like multiple-scattering, reab- 
sorption, or superradiance should be expected. In fact, 
apart from ^ we can also define the analogous of the op- 
tical depth, which for a cubic lattice is x = M^/'^^^. De- 
pending on the values of those dimensionless parameters 
we can predict different phenomena. 

Let us first consider the case of one atom A^ = 1, 
symmetrically distributed through a lattice of M sites. 
We will assume that laser directions are chosen such 
that |ki| = |fco|. The state at time t can be expressed 
as \^{t)) = (l/VM)Ei^i(i)|li,{0}>+Ek^k(t)|0,lk), 
where Ai represents the amplitude at site i. Follow- 
ing similar lines as in the single site example we can 
write Ai{t) — X); r|i-j|^j(i), with rates given by @. 
An analytical solution can be obtained by considering 
a large system (M^'^ ^ 1) such that boundary ef- 
fects are neglected, and periodic boundary conditions 
can be assumed. In that situation, atoms in \a) remain 
in the completely symmetric state with an amplitude 
Aco\\{t) = (l/VM)I]i^i(^) tliat decays with a collec- 
tive rate Fcou — X^n-'^lnl- Then, the system reaches sev- 
eral regimes in which collective effects play an important 
role [18|: (a) If ^ < 1 and x ^ li reabsorption occurs, 
and the decay rate is renormalized to FcoU ^ X^Q- (b) 
If M^/^ > ^ > 1, the dipolar interactions couple near 
neighbors, and x > 1 independent of the size of the sys- 
tem. The rate Fcou scales in the same way as in (a). 
Finally, case (c) corresponds to ^ ^ M^/'^, a situation in 
which every site is connected through dipole interactions 
to all other sites, and Fcoii — MTq. 

We now consider N atoms within M sites. The above 
described regimes still give a valid picture for this situa- 
tion. In particular, we focus on the collective limit (b), 
where ^ > 1. In addition, from here on, we will consider 
laser directions to be such that jkildoAf^/'^ ^ 1. When 
all the atoms are initially in the lattice, we get 



j 

: Er*,_i|(a+afaj)+r|,_j|(a+a?a,), (6) 



dt 




FIG. 3: Rate of atomic emission from a ID lattice with A*' = 
M — 100 atoms. Solid, dashed, dotted and dot-dashed curves 
correspond to ^ = 0.9,1.25,2,3.33 respectively. Inset |17l |: 
Spontaneous symmetry breaking and non-zero steady state 
population. Solid and dashed curves represent z{t) and y{t) 
respectively, evolving according to ([HJ, for N — M = lO'^ 
atoms, A = 0, and uo = 50. 

with rates given by @. Here, a^ = 2a^(Ti — 1, and 
all the operators are evaluated a time t. Let us first ana- 
lyze the atomic emission that occurs for positive detuning 
A > 0. We focus on the rate of emission of atoms in all 
directions, which is given by 7?.(i) « — X]j d{a^)/dt, for 
different values of ^. If sites evolve independently, TZ{t) 
decays exponentially. However, when ^ > 1 and collec- 
tive effects are present, TZ{t) does no longer decay expo- 
nentially and, furthermore, it presents positive slopes at 
initial times. This is shown in Fig. [3] for a ID lattice, 
where it is observed that collective effects occur for ^ > 1. 
This result is obtained with ^ by using the semiclassical 
decoupling (cr/'crfcrj) = (cf )(o'i^crj), that is based on ne- 
glecting atomic quantum fluctuations [9|. Nevertheless, 
the change of sign in the slope can be obtained analyti- 
cally by differentiating Eq. ([6]) at i = without the use 
of any approximation. 

If we now consider negative detuning, A < 0, the rates 
([4]) are purely imaginary, and the system has an effective 
Hamiltonian 



H: 



A<0 



eff 



J2j\i-j\<^l<^i 



(7) 



dt 



where J|i-j| = *r|i-j| is a real and negative quantity, so 
that ([7]) is Hermitian and describes a coherent spin-spin 
interaction of ferromagnetic type. This interaction may 
have interesting applications in the field of quantum sim- 
ulation. Furthermore, for ^ ^ 1 it gives a Coulomb-like 
interaction very difficult to obtain with other techniques. 
We now concentrate in the non-Mar kovian limit, 
where the system also becomes strongly interacting. We 
will consider that all the atoms are initially in the lattice 
(i.e., N — M), and the limit M ^ 1. We use the mean 
field or Hartree approximation [22]. Then, the evolution 
of yit) = Eji^i {t))/M and z(i) = Ej(4(i)>/^ can be 



written as 

dy{t) 
dt 

dz{t) 



dt 



= M / dTG,o\i{t-T)y{T)z{t)- 



dTG,on{t-T)y*{T)y{t) 



= -4MRc 



(8) 



The function Gcoii(^) = X^n C!n{t), with Gn{t) defined in 
([5]). Due to the non-Markovian structure of the equa- 
tions, the mean field approximation here considered pre- 
dicts that the trapped atoms acquire a macroscopic po- 
larization in the steady state. This is shown in Fig. ([3]) 
for A = 0, where we have considered an initial infinites- 
imal polarization, j/(0) = 10~^ (though the steady state 
does not depend on the choice of y(0)), and z(0) = 1. 
The spontaneous polarization of the system, previously 
described in [7| for atoms in PCs, is similar to the sponta- 
neous symmetry breaking described in the semiclassical 
theory of the laser ([22] and references therein). Fig. ([3]) 
also shows that the non-Markovian effects lead to a non- 
zero steady state population, i.e. z*** ^ —1. 

Most of the phenomena described here may be ob- 
served with state of the art experimental setups us- 
ing state-dependent potentials [l| in the Mott insulator 
regime [2| for the lattice atoms and choosing fi <C wq to 
avoid the occupation of other bands. The simplest regime 
corresponds to the a single excitation (A'' = M = 1), 
where decay, non-Markovian effects, as well as the phase 
transition occurring at A = can be observed. Note that 
it is not required to have a single atom in the whole lat- 
tice, as long as the atoms do not interact with each others 
(i.e. ^ <C 1 and initially localized), so that one could eas- 
ily monitor the decay as a function of time by simply 
measuring how many atoms remain in \a). The pres- 
ence of the bound state of the untrapped atoms should 
be visible in the momentum distribution after free ex- 
pansion. By preparing a few atoms each of them co- 
herently distributed among M sites in disjoint regions, 
several copies of the setup consisting in TV = 1 atom 
within M sites could be realized. Hence, we would ob- 
serve collective effects in the decay time for x > 1. For 
an initial Mott insulator state (TV = M) it should also be 
possible to observe superradiant effects by looking at the 
slope of the decay rate for short times, and for A < the 
nearest-neighbor interaction induced via virtual transi- 
tions to the untrapped state, as follows from Eq. ([7]). All 
those phenomena require ^ > 1, i.e. |A| < l/(2?Tido), 
as well as the Markovian limit, FcoU < |A|. Observing 
superradiance for long times (i.e. the whole shape of Fig. 
[3]) may be limited by decoherence effects caused by ran- 
dom magnetic fields which shift \a) and \b) differently. 
Furthermore, to observe coherent interactions in Eq. ([7]) 
beyond nearest neighbors requires ^ 3> 1 , which may also 
be compromised by decoherence effects. A possible way 
around this is to use lighter atomic species, like Li, where 
those conditions are relaxed. On the other hand, the col- 
lecive non-Markovian effects related to the spontaneous 



polarization should be also easy to observe by choosing 
A~0. 

The proposed set-up may be advantageous to observe 
some phenomena with respect to atoms interacting with 
light in a PC. First, it is easily tunable; second, the de- 
tection techniques developed for atoms in optical lattices 
[23| may allow to measure features, like the analogue to 
the photon-atom bound state, that are difficult to mea- 
sure in PCs; third, the optical lattice is a nearly perfect 
periodic potential. In addition, other interesting phe- 
nomena could be explored with the present set-up. For 
example, for x ^ 1 light-matter interface schemes can 
be used to control the emission direction of the atoms, 
or to map the state of trapped atoms into that of un- 
trapped atoms 241 . Besides that, if \b) is affected by a 
trap that is wider than that of \a), other types of inter- 
actions, like the one described by the Jaynes Cummings 
or the so-called Tavis Cummings model, can be imple- 
mented. Finally, the present system may be tuned to 
explore other regimes which have never been considered 
in quantum optics since they could not be reached there, 
like for example the ones in which the initial state of the 
atoms in the lattice is a superffuid, or a Tonks gas. 
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